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^ '. Abstract 

We develop a method to determine the eigenvalues and eigenfunctions of two-boson Hamiltonians 
include a wide class of quantum optical models. The quantum Hamiltonians have been transformed 
in the form of the one variable differential equation and the conditions for its solvability have been 
discussed. Applicability of the method is demonstrated on some simple physical systems. 
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I. INTRODUCTION 



For over ten years there has been a great deal of interest in quantum optical models 
which reveal new physical phenomena described by the Hamiltonians expressed as nonlinear 
functions of Lie algebra generators or boson operators . Such systems have often 

been analyzed by using numerical methods, because the implementation of the Lie algebraic 
techniques to solve those problems is not very efficient and most of the other analytical 
techniques do not yield simple analytical expressions. They require tedious calculations 0, 



However, recently a new algebraic approach, essentially improving both analytical and 
numerical solution of theproblems, has been suggested and developed for some nonlinear 



quantum optical systems 



0, 11, Q. 



Most of such developments are mainly based on 



linear Lie algebras, but it is evident that there is no physical reason for symmetries to be 
only linear. Nonlinear Lie algebra techniques and their relations to the nonlinear quantum 

14, 3,0)- In both cases finite part of spectrum of 



optical systems have been discussed jlH 
the corresponding Hamiltonian can be exactly obtained in closed forms and 
are known as quasi-exactly-solvable(QES) termed by Turbiner and Ushveridze 



;hese systems 



171 ] . Recently 



it has been proven that the single boson Hamiltonians also lead to a QES under some certain 



constraints 



The aim of this paper is to determine quasi-exact-solvability of the two-boson systems 
and discuss their possible applications in physics. As a particular case our model includes 
the solutions of the Karassiov-Klimov Hamiltonian and the Hamiltonian of the systems of 
photons and bosons expressed in a single mode form. These Hamiltonians are _not only 
mathematically interesting but they have also potential interest in physics 



20 



21 



22 . 



The paper is organized as follows: In section 2, general form of the Hamiltonian has 
been constructed and its solution by using the invariance of the number operator have 
been discussed. This section includes some physical examples. In section 3 we present a 
transformation procedure in order to solve a wide range of the Hamiltonian. In section 4 
the Hamiltonian is transformed in the form of the one variable differential equation and the 
conditions to obtain its eigenvalues and eigenfunctions are discussed. Finally, in section 
5 we comment on the validity of our method and suggest the possible extensions of the 
problem. 
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II. TRANSFORMATION OF THE TWO-MODE BOSONIC SYSTEMS IN THE 
FORM OF SINGLE- VARIABLE DIFFERENTIAL EQUATION 

Two mode bosonic Hamiltonians play an important role in nonlinear quantum optical 
systems. The Hamiltonians of such systems can be generalized as follows: 

H = ^« mi , m2 , m3 , m4 (a+) mi («i) m2 (4) m3 («2) m4 (1) 

mi 

where q.{ is a constant and rrii determines order of the interaction. The boson creation, a\, 
a 2 , and annihilation af, operators obey the usual commutation relations 

[ai,a 2 ] = [ai,aj] = [a 2 ,a|] = [a|,aj] = 0, [ai,af] = [a^aj] = 1- ( 2 ) 

We assume that the conserved quantity of the physical system described by the Hamiltonian 
H, correspond to the operator: 

K = safcti + pa 2 a 2 . (3) 

Clearly the Hamiltonian conserves the number of particles in the system, when H commutes 
with the K. Meanwhile, we have mention here, in general, conserved quantity of (j2J) type 
appears more general contexts than the number operator. Indeed a convenient quantum 
mechanical treatment of the m th order interaction of the two-boson model exploits the 
SU(2) algebra structure of the problem, by expressing the Hamiltonian as a polynomial of 
the angular momentum operators whose degree depends on the degree of interaction. In 
this case the conserved quantity is the total angular momentum, J = j(j + 1), which can be 
related to the number operator j = N/2, where N = a^ai + a\a<i. Therefore the conserved 
quantity J crucially depends on the conservation of the number of particles. The classical 
motion of the particle takes place in the space of angular momentum on the sphere of radius 
j(j + 1). To have a physical insight in to K, from other perspective, when the motion of the 
particle takes the place on an ellipsoid the conserved quantity K takes the form (jSJ). The 
operator K and bosonic operators satisfy the commutation relations 

[K, af] = saf, [K,ai] = —sai, [K, ] = pa^", [K, a 2 ] = — pa 2 . (4) 

If K and H commute, then they have the same eigenfunction. Therefore it is worth to 
seek the conditions for the commutation of H and K. This can easily be done by using the 
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commutation relations (J2J) and (@J) and we obtain the following relation 

[K,H] = ^[(s(mi-m 2 ) + p(m 3 - m 4 )]« mim2m3m4 (a+) mi (a 1 ) m2 (a+) m3 (a 2 ) m4 . (5) 

mi 

It is obvious that the constant of motion K and H commute when the following set of 
equation is satisfied 

s(mi — m 2 ) + p(m 3 — m 4 ) = 0. (6) 
The action of the operator K on the state \ni,n 2 ) is given by 

K\n 1 ,n 2 ) = (sn 1 +pn 2 )\n 1 ,n 2 ) . (7) 

In this section we seek the solution of the eigenvalue equation (JJJ) in the Bargmann-Fock 
space. The usual realization of the bosonic operators in Hilbert space is given by 

4 = ai = 7f(^ + Xl ) (8a) 

lid \ lid 

1- To 1 ' fflo = 



7T 



" 2 = tA~s^ +X2 )' tt2 = 7lU^ + I V (8b) 

The operators can be transformed to the Bargmann-Fock space by introducing the operator 

exp - (a\ + a\ 2 + a\ + aj 2 ) (9) 
L 8 

and the similarity transformation 

1 1 
bi = TaiT -1 = —= (a~i — ai) , = Ta^r -1 = —= (af + ai) 

b 2 = ra 2 r _1 = -^= (a 2 - a 2 ) , b 2 = Ta^r" 1 = — (a 2 + a 2 ) . (10) 
v2 v 2 

The Bargmann-Fock space differential realizations of the operators take the form 

&! = -—; bj = xi] 6i = — ; 6j = x 2 (11) 
axi ax 2 

Thus, in the Bargmann-Fock space the eigenvalue problem (J7J) leads to the following solution: 

1>{x ll x 2 )=a%<j>(x 2 xSy (12) 

where k is given by 

k = n 1 + ^n 2 . (13) 
The eigenfunction of the Hamiltonian can be obtained from the relation 

|ni,n 2 } = Y~ x ^{xx,x 2 ). (14) 
In the following section the application of our procedure will be discussed explicitly. 



A. Hamiltonian of the second Harmonic generation and generalization of the solv- 
ability conditions 

The Hamiltonian (JTJ) includes various physical Hamiltonians, second and third harmonic 
generation effective Hamiltonians. This gives us an opportunity to test our approach because 
those Hamiltonians have been studied in literature. Consider the following Hamiltonian: 

H = uJiafai + tu 2 a 2 a 2 + /t(a^) 2 a 2 + ~Ka 2 (ai) 2 (15) 

which can be related to the Hamiltonian ([TJ, when the parameters 

C^mim2,r7i3,rri4 (16) 

except that 

"1,1,0,0 = U X , a 0)0 ,l,l = ^2, "2,0,0,1 = K, "0,2,1,0 = « (17) 

The condition © is satisfied when p = 2s. The constants u\ and uo 2 are the angular fre- 
quencies of two independent harmonic oscillators characterized by annihilation and creation 
operators, k and ~k are the coupling coefficients that determine the strength of the interaction 
of the oscillators. In the following we use the Bargmann-Fock realization, where creation and 
annihilation operators a+ and a.; are replaced by bf and fej, respectively. The eigenfunction 
()12|) is of the form 

ip(x 1: x 2 ) = (xi)VW (18) 

where z = x 2 x l 1//2 . The solution of this system describes a quantum mechanical state of 
H provided that <f>(z) belongs to the Bargmann-Fock space. The eigenvalue equation of the 
second harmonic generation Hamiltonian can be written as 

H^(x u x 2 ) = E^( Xl ,x 2 ) (19) 

Insertion of (JTHj) into (fT9"j) yields the following differential equation 



dz 2 



+ (k + z(uj 2 — 2lo\ + 2kz(3 — 2fc))— — \- uj 2 + kui + Kzk(k 

( I .-_ 



1) — E 



0. 

(20) 
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According to Turbiner|l7j] the Hamiltonian (j20j) is quasi-exactly solvable. The eigenvalue 
equation (j2*Uj) can be obtained as follows. Let us assume that the function <p(z) is a polyno- 
mial in z with the coefficients being functions of energy: 

oo 

ct>{z) = Y,Pm{E)z m . (21) 

m=0 

When inserted <j>{z) in (J2U|) then we obtain the following three-term recurrence relation: 



K(k-2m)(k-2m-l)P m+1 (E)+ (22) 
(u 2 + kuji + m{uj2 — 2u{) — E)P m (E) + KmP m -i(E) = 0. 

The function P m (E) terminates when m > k/2, and therefore the roots of the P m (E) belong 
to the spectrum of the Hamiltonian. For the corresponding problem, Beckers, Brihaye 
and DeberghQ] have obtained three different Hamiltonians in the framework of nonlinear 
algebra. Two of the Hamiltonians are PT- symmetric and they related to the sextic harmonic 
oscillator. The third Hamiltonian is exactly the same one in (j2"Uj) which was obtained in 
a different method. The Hamiltonian (|20j) can also be transformed in the form of the 
Schrodinger equation by changing variable z = —l/~K,y 2 and defining the wavefunction 

<f)(z) = e-^ w{y)dy ^(y) (23) 

where 

W{y) = * + 0* -**)V _ ^! ( 2 4) 
y 4 4 

Substituting (|2~3j) in (|2"U|) we obtain the following equation: 



H = + - ({2k + 5)w 2 - 2o;0) + 4 ((^ - 2^) 2 - 4^(2A; + 3)) y 2 - 

2 ay z 4 lb 

\kk(u 2 - 2u x )y 4 + ^ 2 ^V- (25) 
o lo 

The last equation is the Hamiltonian of the sextic Harmonic oscillator. Consequently the 
procedure given here can be applied to obtain eigenfunction and eigenvalues of the var- 
ious physical Hamiltonians. The validity of the procedure depends on the choice of the 
Omim2,m 3 ,m 4 -One can easily be obtain various physical Hamiltonians by appropriate choice 



b 



of ctm 1 m2,m 3 ,m 4 and by considering the conditions given in (JSJ). For instance one can easily 
obtain n th harmonic generation Hamiltonian by setting: 

a l, 1,0,0 — k'lj a 0,0,l,l — a n,0,0,l — K ! a 0,n,l,0 = K (26) 

and the remaining values of the a mim2im3]m4 = 0. Then the Hamiltonian (0) takes the form 

H = Uiafai + u 2 4 a 2 + n(af) n a 2 + K4( a i) n - (27) 

In the following section we continue to explore the solution of the Hamiltonian (JI)), in the 
Bargmann-Fock space by developing a similarity transformation procedure. 

III. TRANSFORMATION OF THE BOSONIC OPERATORS 

In this section we present two transformation procedure to obtain the conditions of the 
solvability of the Hamiltonian (JI}. The following transformations allow us to study a wide 
range of physical systems. Let us introduce the following similarity transformation induced 
by the operator 

S = (a+)"^ ai (28) 

where rj is a constant and it can be determined by considering the transformation of the 
number operator (J5J). The operator S acts on the state \rii, n 2 ) as follows, 

S |ni,n 2 ) = (4) vni \n 1 ,n 2 ) = \ — — — rr \n 1 ,n 2 + v n i) ■ ( 29 ) 

V \ n 2 + V n iY- 

Since and a 2 commute, the transformation of a\ and af under S can be obtained by 
writing 4 = e b , with [ax, b] = [4, b] = 0, 

Sa+S- 1 = e ^ata la + e -vba+ ai = +( a +)»? 

SaxS- 1 = e nbaia laie - V ba+ ai = a ^ a +yv ( 3 Q) 

and transformation of a 2 and a 2 as follows 



sats- 1 = (4)^^4(4)-^^ = 4 

Sa^- 1 = (4Y a > 1 a 2 (4Y^ ai =a 2 - r]4a 1 {a+)-\ (31) 
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In a similar manner, we prepare the other similarity transformation which is useful to study 
QES of the Hamiltonian (0), by introducing the following operator: 



T = aT 1 ai (32) 
The operator T acts on the two-boson state as 



T\n 1 ,n 2 ) = a 2 ni \n 1 ,n 2 ) = x l -, — — -r \n 1 ,n 2 - cmi) . (33) 

[n 2 — anijl 

Since a\ and a 2 commute, the transformation of a± and af under S can be obtained by 
letting a 2 = e c with [a 1; c] = [af, c] = 0, 



Ta+T~ l = e ^i"i a + e - Qca i a i = a+(a 2 ) Q 

TaiT" 1 = e ^a+«i aie -^a+ ai = ai ( fl2 )-« (34) 



The transformation of a 2 and at is as follows: 



J -/ = ^2 a 2 a 2 = a 2 "r cnajaia 2 ■ 

1 a 2 l = a 2 a 2 a 2 = a 2 (-35) 

These two transformations lead to the conditions of the quasi-exact-solvability of the Hamil- 
tonian (fT|). 

IV. SOLVABILITY OF THE HAMILTONIAN 

In this section we discuss the solvability of the Hamiltonian (0). The conserved quan- 
tity K of the physical system describes the states of the corresponding Hamiltonian. The 
transformation of K under the operator S is given by 

K' = SKS~ X = (s — pr))afai + pafa 2 (36) 

The Hamiltonian $IJ is characterized by the total number of a\ and a 2 bosons when the 
conserved quantity of the physical system K commutes with the whole Hamiltonian. In the 
transformed case it is only the number of a 2 bosons that characterize the system under the 
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condition rj = s/p. According to (|36|) the representation is characterized by a fixed number 
at 0-2 = k. Therefore in the transformed form the Hamiltonian can be expressed as one boson 
operator a%, when the condition JHJ) is taken into consideration. The transformed form of 
the Hamiltonian (JTJ can be written as: 

H = SHS- 1 = ^« mi , m2 , m4 (a+) mi («i) m2 (^ - -afa^ (37) 

mi 

The difference between JIJ and (J3~T|) is that while in the first the total number of a\ and 
a 2 bosons characterize the the system, in the later it is only the number of a 2 bosons that 
characterize the system. Therefore the representation is characterized by a fixed number k 
and in (|37|). the Hamiltonian is expressed in terms of one boson operator a\. The transformed 
Hamiltonian H, in the Bargmann-Fock space, which plays an important role in the quasi- 
exact solution of the equation (0). It can be transformed in the form of the one dimensional 
differential equations in the Bargmann-Fock space when the boson operators are realized as 

ai = — , at = x. (38) 
ax 

The basis function of the primed generators of the system is the degree of polynomial of 
order k, 

P n (x) = (x\x\--- } x k ). (39) 
Action of (|37|) on the (}3*9"j) . in the Bargmann-Fock space can be written as: 

HP n (x) = J2 P n(E)x n (40) 

where E is the eigenfunction of the Hamiltonian and the polynomial P(E) can be written 

as 

V a mi m2 m4 (k - -n) m \ nl P n+mi _ m2 (E) - EP n (E) = (41) 

L — 4 ' ' p [n — mo)! 

mi v 1 

The wavefunction is itself the generating function of the energy polynomials. The eigenvalues 
are then produced by the roots of such polynomials. If the Ej is a root of the polynomial 
Pj + i(E), the series (f4*Tj) terminates at j > k- and Ej belongs to the spectrum of the 
corresponding Hamiltonian. The eigenvalues are then obtained by finding the roots of such 
polynomials. 
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The constant of motion K characterize the system can be transformed, by the operator 
T is given by 

K' = TXT' 1 = {s+ prj)afa 1 + pa+a 2 . (42) 

The Hamiltonian ((TJ can be characterized, in the transformed case, 77 = —s/p. Thus accord- 
ing to (}4*2*|) the representation is characterized by a fixed number a\a 2 = k. Therefore the 
transformed Hamiltonian includes one boson operator a 1; when the condition (JUJ) is taken 
into consideration, and it can be written as: 

H> = THT- 1 = Va mi m2 m M) mi {aiT 2 {k + ^'"^ + -a+a^ (43) 
/-^ ' ' pm 3 p 

The Hamiltonian can be expressed as one dimensional differential equation in the Bargmann- 

Fock space. In order to solve the Hamiltonian H' we can follow the procedure given (J38|) 

through (jHJ). In this case the polynomial function is terminated when k is constrained to 

sfmi — mo) s . . 

k > — y — ^ - -n 44 

pm 3 p 

Consequently we have obtained two classes of Hamiltonians whose spectrum can be obtained 
(quasi)exactly. 



V. CONCLUSION 



In this paper we have prepared a general method to obtain the solution of two boson 
Hamiltonian. By using either solution of number operator or similarity transformation, we 
have been able to provide a QES of the various physical Hamiltonians. Furthermore, it 
has been given that two boson Hamiltonian can be reduced to single variable differential 
equation in the Bargmann-Fock space. It is also important to mention here that the methods 
given here can be used to solve higher order differential equations. 

The method given here, can easily be extended to solve the Hamiltonians that include 
multi-boson or fermion-boson systems. This extension leads to the solution of the various 
physical problems, such as Pauli equation, Rabi Hamiltonian, Jaynes-Cummings Hamilto- 
nian. 
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